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Abstract 

We carry out a detailed analysis of the scalar sector of the nonminimal supersymmetric standard 
model with lepton number violation, and study the constraints imposed on it by the stability of the 
electroweak symmetry breaking vacuum. The model contains a trilinear lepton mumber violating 
term in the superpotential together with the associated supersymmetry breaking interactions 
which can give rise to neutrino masses. We evaluate the mass matrices for the various boson and 
fermion modes and then discuss the effect that the lepton number violating interactions have on 
the mass spectra using a phenomenological precription to implement the experimental constraints 
on the light neutrinos mass matrix. We also discuss qualitatively the conditions on the lepton 
number violating parameters set by the unbounded from below directions, and from the absence 
of the charge and color breaking minima in this model. 
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I. INTRODUCTION 



There are suggestive hints of lepton number violation in the observation of neutrino 
flavor oscillations [1] . The link between the new physics beyond the Standard Model (SM) 
and lepton number violating interactions has so far been realized in two main theoretical 
frameworks. The first option, the see-saw mechanism [2], uses physics at the grand unified 
scale, with lepton number violation represented by the non-renormalizablc interaction of 
dimension 5, >Ce// = \{LHy){LHy), where L and are the lepton and Higgs doublets, 
respectively. The second option involves minimal supersymmetric standard model (MSSM) 
with broken (generalized) lepton parity, and uses bilinear and/or trilinear lepton number 
violating couplings in the superpotential, with lepton number violation occurring at the 
weak scale [3]. 

Recently, we have pointed out [4] that in the context of low energy supersymmetry there 
is an attractive option of generating neutrino masses in the nonminimal supersymmetric 
standard model (NMSSM) with trilinear lepton number violation. Unlike the minimal su- 
persymmetric standard model (MSSM), the superpotential of the NMSSM is scale invariant. 
What is perhaps even more interesting is the presence of an additional lepton number vio- 
lating trilinear coupling in the superpotential of NMSSM which has no analog in the MSSM. 
In the NMSSM with lepton number violation we have a situation which is different from the 
one that arises in the see-saw mechanism or the bilinear lepton number violation in MSSM 
in that no dimensionful mass parameters (large or small) are introduced. 

There are several reasons for studying the nonminimal supersymmetric standard model. 
First, and foremost, it is the simplest supersymmetric extension of the Standard Model in 
which the electroweak scale originates from the supersymmetry breaking scale only. Second, 
and as mentioned above, this model breaks lepton number and i?-parity explicitly through 
a unique trilinear superpotential term [5, 6]. Third, the model can successfully explain the 
pattern of neutrino masses with normal type of hierarchy [4] without invoking either a large 
or a small mass parameter. And last, but not the least, its enlarged Higgs sector may help in 
relaxing the fine-tuning and little hierarchy problems of the MSSM [7] , thereby opening new 
perspectives for the Higgs boson searches at the next generation of high energy colliders. 

In addition to the familiar Yukawa interaction superpotential for the three generations 
of quark and lepton superfields, and for the Higgs bosons superfields, Hu, Ha, the model 
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couples the set of down type Higgs boson and sleptons, La — {Ha, Li) {A — d,i; i — 1,2,3), 
with the Higgs singlet superfield S through the superpotential W — \aLaHuS — ^S^, where 
\a = (A, Ai) and k are dimensionless parameters. The familiar soft supersymmetry 
breaking trilinear scalar couplings and mass terms introduce the dimensional parameters 
^4^^ = {Ax,Aj^,), Afi, and fn^. Whereas the fundamental parameters in the action 

are A^, ^Aa' ''^\ab^ vacuum solution is specified by the dynamical parameters corre- 
sponding to the vacuum expectation values (VEVs) of the electrically neutral scalar fields 
at the electroweak scale, Vu —< >, va — {vd —< >,Vi —< i>i >), which include the 
sneutrinos together with the neutral components of the Higgs bosons. We use a covariant 
notation for the down-type Higgs boson and slcpton supcrficlds, and their couplings, in order 
to illustrate the invariance of obscrvables under the SU (4) group of transformations acting 
on the 4D column vectors La and A^, ^Aa' matrix m| Strong correlations usu- 

ally exist between the fundamental and the dynamical parameters. Thus, the experimental 
observations for the light neutrino mass matrix tightly constrain certain combinations of the 
parameters which characterize the inequivalent field directions distinguishing the down- type 
Higgs boson from the sleptons, and hence the interactions which conserve lepton number 
from those which violate it. One may hope that imposing the necessary restrictions on the 
misalignment type parameters would still leave some freedom in the parameter space so as 
to allow the lepton number violating interactions to significantly affect the physics of the 
NMSSM. There is a close analogy here with the situation that obtains in the MSSM with 
R-parity violation [8-15] . 

At the tree level, the size of the neutrino Majorana mass matrix is knownn to be controlled 
by deviations from the A^ oc va alignment, such that the predicted single no n- vanishing neu- 
trino mass is proportional to J2i A^^, where the effective misalignment parameters are defined 
by the linear combinations X'- = X^ — Xvi/va- We have previously shown [4] that the suitably 
constrained NMSSM predicts an acceptable neutrino mass matrix of a normal hierarchy 
type so that the tentative hmits on the heaviest neutrino mass, m^^ < (10 MeV — 1 eV), 
set bounds on the effective parameters given by X[ < (3. 10~^ — 10~^). Bounds of similar 
type were found for the couplings Aj in the basis Vi = in an attempt to fit an extended 
version of the model to the neutrino mass differences and mixing angles [16]. It has also 
been observed that ratios as large as Aj/A ~ 10~^ might be allowed at the price of fine 
tuning the NMSSM parameters [17]. As far as one- loop contributions to the neutrino mass 
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matrix is concerned, the situation is more complex due to the larger number of unknown pa- 
rameters and the presence of several sources of contributions arising in perturbation theory. 
Thus, to supress the contributions from the one-loop amplitude involving the exchange of 
scalars and gauginos, it is necessary to impose strong constraints on the effective parameters, 
B'- = Bi — Bvi/vd [Ba = {B, Bi) = (A^^.Aj, ^A''^)] controlled by the Ba oc va alignment. 

Supersymmetric models, including the nonminimal supersymmetric standard model, con- 
tain a large number of scalar fields, such as the scalar partners of quarks and leptons, which 
are not present in the Standard Model. As a result, the scalar potential of supersymmetric 
models is fairly complicated. Such a potential can have minima which could lead to the 
breakdown of color or electric charge conservation through the vacuum expectation values 
of squark and slepton fields. The condition that the global minimum of the theory does 
not violate color or electromagnetic gauge invariance provides useful constraints on the pa- 
rameter space of the underlying supersymmetric model, as has been well illustrated in the 
context of the minimal supersymmetric standard model [18-22] and its extended version 
with broken R parity symmetry [23, 24]. Significant bounds can be obtained not only on 
the value of the trilinear soft supersymmetry breaking parameters, but also on the values 
of the bilinear soft scalar supersymmetry breaking parameter, as well as scalar and gaugino 
masses, respectively. 

The possibility that one can reproduce the lepton number violating observables without 
imposing strong restrictions on the coupling constants Aj motivates us to examine how 
the scalar sector of the NMSSM is affected by treating Aj as free parameters. One may 
expect interesting implications regarding the stability of the regular vacuum solution and 
the properties of the neutral and charged scalars. While performing the calculations within 
a field basis independent formalism [8-10, 12-15] is highly desirable, this is very tedious and 
is not particularly illuminating. It requires an unusually large amount of effort which is not 
warranted at the present stage of our study. We, thus, choose to pursue the calculations 
within a basis dependent formalism. In principle, there is a choice of basis for the neutrino 
superfields in which the sneutrino VEVs can be set to zero at the expense of transferring 
through a superfield redefinition the couplings LiH^S to the trilinear lepton number violating 
couplings of the matter superfields. Although the calculations are greatly simplified in the 
basis Vi — 0, there is a risk that this choice biases the parameter space exploration, as has 
clearly been pointed out in the above discussion regarding the neutrino mass matrix. In 
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fact, in the basis choice Vi — 0, the couphngs Aj are so strongly suppressed that the issue of 
estabhshing any observable effect on the scalar sector would be closed from the very start. 

In the present work, we shall adopt the view that the lepton number violation phe- 
nomenology allows the couplings Aj to cover a wide range of variation, < |Aj/A| < 1. 
Thus, we treat all the A^ on the same footing, while requiring the soft and dynamical pa- 
rameters to satisfy the near alignment properties, va oc Ba and va oc A^, in agreement 
with the observations. Following a phenomenological approach, we treat Aj as free param- 
eters, and determine Vi and Bi by restricting the linear combinations X[ — Xi — Xvi/vd and 
B'j^ = Bi — Bvi/vd- Since v^, Vi are determined by minimizing the scalar potential, this pro- 
cedure implies a fine tuning. However, one can use symmetries to impose in a natural way 
a preferred direction in the La field space along which the various parameters get aligned. 
While the horizontal flavor discrete symmetries [25, 26] indeed have the ability to enforce 
an approximate A^ oc va alignment, this property is robust only in the case of hierarchically 
small ratios, ^ ~ which is not exactly what we want. By contrast, the supergravity 
inspired models of the type discussed by Nilles and Polonsky [27], in which the XaLaH^S 
couplings arise from non-renormalizable interactions controlled by a spontaneously broken 
U{1)r type symmetry, are compatible with all the A^ having the same size. The alignment 
is realized dynamically via universal type boundary conditions on the soft parameters, and 
remains typically stable under the renormalization group [27, 28]. Our main purpose is to 
study the impact of the lepton number violation on the electroweak symmetry breaking vac- 
uum of the nonminimal supersymmetric standard model consistent with the experimental 
measurements for the light neutrinos. In the first stage, we discuss its stability by testing 
for the occurrence of tachyons (negative eigenvalues of the squared mass matrices) for the 
charged and neutral scalar modes. In the second stage, we analyze the dangerous directions 
in the scalar field space along which the scalar potential is unbounded from below (UFB) 
or develops deeper charge and/or color breaking (CCB) minima. Distinct methods must 
be used at these two stages, which both involve the non-supersymmetric couplings in an 
essential way. 

This paper is organized as follows. In Section II we describe the nonminimal supersym- 
metric standard model with lepton number violation, and set up our notation and con- 
ventions. Complementary definitions and results are summarized in the Appendix A. In 
Section III we discuss the conditions on the superpotential couplings and supersymmetry 
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breaking couplings required to exclude the tachyon scalar modes for the regular vacuum. 
In Section IV we explore the conditions on the superpotential and supersymemtry breaking 
couplings implied by the lifting of UFB directions, and the removal of CCB minima. In 
Section V we summarize our main conclusions. 

II. THE NMSSM WITH BARYON AND LEPTON NUMBER VIOLATION 

In this section we summarize the basic features of the superpotential and the associated 
soft breaking terms for the nonminimal supersymmetric standard model with baryon and 
lepton number violation, and establish our notation and conventions. Some of the details 
are further discussed in the Appendix A. 

A. Superpotential and soft breaking terms 

The superpotential of NMSSM with baryon and lepton number violation is characterized 
by the scale invariant supepotential 

W - Wrpc + Wrpv, (II. 1) 

where Wrpc and Wrpv arc the baryon and lepton number (and i?-parity) conserving, 
and baryon and lepton number (and i?-parity) violating contributions, respectively to the 
superpotential [4]. These contributions are written down explicitly in Appendix A. In 
the absence of lepton number conservation, there is no distinction between the down-type 
Higgs {Hd) and the lepton superfields (Lj) as they transform identically under the SM gauge 
group. We can, therefore, employ a four vector notation for the down type Higgs and the 
lepton superfields, and the Yukawa couplings 

La = {Hd.Li), 

Xa = (A, A,), [A = d, i- i = l,2,3], (II.2) 

which is given in detail in Eq.(A.4) of Appendix A. The quark and lepton generation 
indices are denoted by the letters i, j, • • •, and the trilinear couplings obey the antisymemtry 
property, A^Bfe = — Ab^A;- In this notation we can write the complete superpotential of 
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NMSSM with baryon and lepton number violation as 

W = X;,H^Q,Ul + \\ABkLALBEl + >^AjkLAQ,Dl + \>^l^,UtD'^pi 

+ \aLaHuS - "^S^ (II.3) 

The supersymmetric contribution to the Lagrangian of NMSSM with baryon and lepton 
number violation generated by the superpotential (II. 3) can be obtained by a standard pro- 
cedure. To this we must add the supersymmetric contribution from the D-iervas generated 
by the gauge interactions. In addition to the supersymmetric part, the Lagrangian consists 
of soft supersymmetry breaking terms, which include soft trilincar scalar couplings and soft 
masses for all scalars and gauginos, respectively. The part of the soft supersymmetry break- 
ing scalar potential of NMSSM containing trihnear scalar couplings can be written as (note 
that it is — y that occurs in the Lagrangian) 

l^^°f^(trilinear) = V|^^^(trilinear) + V|^i^(trilinear), 



— --^^ABk^ABkLALsEl - A^^^Xj^j^LAQjDl - A'^^y^^HuQjU^ 
-l^kKjkU^Dpl - A-,^~XaLaH^S - + H. c, (II.4) 



and the part containing the bilinear mass terms for the scalars and gauginos can be written 

as 

T/^oft(mass) = (mass) + (mass), (II.5) 

with the lepton number conserving RPC terms given by 

■S0ft/_„„„\ _ 2 i/) |2 I 2 \fTC\2 I 2 \f\c\2 , 2 | r |2 , 2 | ™|2 



v^^^^(mass) = m^.|gi|' + m'^e|t/rr + ^^clAT + ^i.l^ir + ^il^, 
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- ( \MaKK+ H. c), (n.6) 

a=3,2,l ^ 

and the lepton number violating terms given by 

V|oft(mass) = ml^^HlL, + i^l^H^HaH^ + i^l^^LiH^+ H.c. (II.7) 

We have used the convention in which the repeated indices are implicitly summed over and 
have suppressed the SU{3) color quantum numbers of the quark and squark fields by set- 
ting, for instance, Qf = {U",Df) — > = {Ui,Di). Although we have represented the 



fields mass mixing by the general terms, m-^.^-LlLj and ^^V^Lb, for the sleptons and 
for the down Higgs boson and sleptons, we neglect thoughout the present work the intergen- 
erational mixing of sleptons by assuming that are mass basis fields with tti^^j = Trij^ Sij 
while retaining the off-diagonal mass parameters, fn^^ i_ = ni^di 0- Since supersymmetry 
breaking respects the SM gauge symmetry, we must assign the same soft mass parameters 
for the two members of electroweak doublet fields, as is explicitly done for the expression 
of the slepton mass terms, viz. m| iLjP = m~ (li/jP + lej^). The choice of the sign of M2 
in the chargino mass matrix conforms with the convention most often used in the litera- 
ture. We have also included above, for generality, the holomorphic mass mixing terms with 
the parameters fi\ = {fi'jj tt = fi'i , = /i? ), although the soft bilinear operators 

H^^LaHu are, in principle, disallowed, since the same discrete symmetry that forbids the 
bilinear superpotential terms /iaLaHu should also forbid the bilinear mass mixing operators. 
Mass mixing terms of same structure would still arise as effective contributions through the 
soft trilinear operators, —A^^Xa < S > LaHu, with //^^ = — ^4^^^^ < S > . 

The scalar potential receives contributions from the super symmetric as well as the soft 
supersymmetry breaking interactions involving the electrically neutral and charged complex 
scalar fields. The supersymmetric contribution arises from the superpotential (II. 3) (the 
F-terms), and from the I?— terms. The Higgs boson and slepton F-term contributions to 
the scalar potential can be written as 

yp(Higgs, sleptons) = {Wnof + \Wh+\^ + \Ws\^ + + \WeJ^ + \Wec\\ (II.8) 

where 

Who = -X^^Ujul + Xavax, Wfj+ = X^^djul - Xacax, (II.9) 
Ws = XaVaVu - XAeAV+ - kx'^, (11.10) 
= XABkeBel + X'^j^djdl + XAVuX, (H-H) 
We^ = -X'Aji^Ujdl - XsAkVBel - Xav+x, (11.12) 

We'r = XABiVAes, (11.13) 

with various vacuum expectation values denoted by 

VA = (Vd =< >, Vi =< Ui >), CA^ (V_ ^< >, Ci =< Ci >), (11.14) 



< ~e1 >, =< H+ >, =< Hi > . (11.15) 



The F-term contribution of the squarks to the scalar potential can be written as 



l^F(squarks) = jW^^J^ ^ |^^j2 ^ |w |2 ^ |w ^|2^ 



(11.16) 



where 



with the vacuum expectation values of the squarks defined by 



Ui ^ <Ui >, di =< A >, < =< U- >, =< £>■ > 



The D-term contributions to the scalar potential can be written as 



(11.17) 
(11.18) 



(11.19) 



(11.20) 



with 



-^kf +^Kr-kr-ie.r+2Kr)', 

„2 r 

TrS'C/(2) y2 /i |2 I |2 I I |2 I |2 I I |2 i |2 , i 



★ J |2 



+4|t;>_ + t;*ei + i;„ + u*d, 



V, 



SUi3) _ 93 



D 



r-Kr-i^?r 



The complete soft supersymmetry breaking scalar potential is given by 

y-/* = ^|^(|„^|2 + |^^|2)+^2.^^|^c|2^^|^^|^c|2^^2^(|^^|2^|g^|^ 



+ 
+ 



(11.21) 

(11.22) 
(11.23) 



-^x^^a{vaVu - eAV+)x - -A^nx^ + A'-^W^^iujV^ - djV+)ul - A^^^XABkVAeBel 



-^%k>''AjkivAdj - eAUj)dl - ^A^-t^Ml^t^uld'^dl + H. c. 



(11.24) 



We recall our convention of summing over the repeated indices for the (suppressed) color 
indices and for the squark and slepton generation labels, j, /c, • • • = 1, 2, 3. Assuming that 
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the squark fields point in a fixed direction in the color space, no summations over the color 
indices would be present in the above formulae. Furthermore, assuming that the VEVs for 
the squarks and sleptons of different generations affect the scalar potential independently of 
each other, allows one to disregard the summations over the generation labels. 

In order to discuss NMSSM in a somewhat realistic manner, it is necessary to include the 
one-loop contributions to the scalar potential. We shall make the usual assumption [29, 30] 
that the scalar potential is dominated by the top quark and squark modes whose explicit 
contributions are given by 



V, 



loop 



327r2 



(11.25) 



where mt — XtVu is the top quark mass and m|., {i — 1,2) denote the squared masses of the 
top squarks. We do not discuss this point any further, as detailed discussions can be found 
in Refs. [29-32]. 



B. Symmetry Constraints on the Parameter Space and Choice of Free Parameters 

The scalar sector involves the scalar field components of the doublet and singlet Higgs bo- 
son superfields, and the lepton and quark superfields. The parameter space of the NMSSM 
with lepton number violation consist of the gauge couplings Qa {a — 3, 2, 1), the supersym- 
metric couplings A, k, A"j., A^;;., X^^, Xi, Xijk, X[jf., the soft supersymmetry breaking cou- 
phngs A-^^Xa, A^k, A\j^i,XABk, A^j^X'Ajk, ^"fe^"*;, and the soft supersymmetry breaking 
mass parameters, m% , mf , rrif^^, m^j , mjj , m? ,m|, m^^, f , fj,\ . We summarize 
below some useful definitions and abbreviations used in this paper: 



m 



w 



Hill, ^1\ 1 HI ' HI t 1 , 2 I 2\ 



2 V w ' a/ ' zj 2 
2 I 2 

1. Q ri2 9\-^ 92 2 2,-2 -2 2,2 /tt 

tan/?=— , — - — , V =v^^v^, v^ = v^ + v^. (11.26) 

The clcctrowcak symmetry breaking scale has a numerical value v = ~ 174 GeV. 

Without loss of generality, the parameterization of the NMSSM can be simplified by using 
the symmetries of the action and the independence of observables under phase redefinitions 
of the fields. The invariance of the Lagrangian under the SU{2)l x U{1)y gauge symmetry 
allows eliminating four real field degrees of freedom, independently of the structure of the 
Yukawa couplings. We choose the convention where the VEVs of the up-type Higgs boson 
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electroweak doublet are set as =< >= and =< >£ R^- With the choice 
Vj^ — 0, the minimization with respect to the field becomes trivial. Next, using the 
scalar potential invariance under phase redefinitions of the fields H^, U'^,Q and the pair 
of fields Hj, D'^, one can make the following choice [33] for various parameters 

kA^ e R+, XAx e A"A„ e Ui e di e (11.27) 

For completeness, we also observe that in the presence of lepton number violation the 
NMSSM still satisfies a Peccei-Quinn symmetry in the hmit k ^ and a R type symmetry 
in the limit Ax — > 0, A^^, 0, where the U{1) pq and U {1)r groups are defined by the assign- 
ment of charges for {Q, U\ D% L, E\ H^, H^, S), Qpq = (-1, 0, 2, -1, 2, -1, 1, 0), Qr = 
I (3, 0, 2, 1, 4, 1, 3, 2). In the limit x — > oo, with Xax — —/ia, i^x, A^^Xax — —BAfJ'A fixed, 
the physical observables must reduce to those of the MSSM with bihnear it!-parity viola- 
tion [24]. 

In order to obtain significant contributions to the scalar sector observables from the lep- 
ton number violating interactions, some subset of the parameters Aj, A^^^, rn?^^j^ , Vi should 
assume large enough values. Whether this can be achieved while imposing at the same time 
highly suppressed contributions to the neutrino mass matrix is closely related to the formal 
symmetry under the SU (4) group symmetry in the field space La — {Hd, Li) . The freedom 
with respect to the choice of La field basis entails that physical observables can only de- 
pend on the SU{A) invariant combinations of the parameters, A^, Ba — A^^Xa, m| va, 
transforming as vectors or tensors. The physical observables can only depend on the sin- 
glet combinations. At the quantum level, the field basis independence holds only after 
summing all contributions at a given order of the perturbation theory. The basis inde- 
pendent formalism can be developed along similar lines as for the MSSM with R-parity 
violation [8-10]. Since Hd is distinguished from Lj by not having lepton number, it fol- 
lows that the lepton number violating observables can only depend on the invariants of 
the angle type characterizing the inequivalent directions assigned to Ha- Thus, the tree 
level contributions to the neutrino mass matrix involve only the misalignment parameter 
|A A — ^J2A,Bi^ ^ ''^)'ab = \ J2a,b{^aVb — XbVaY, while the one-loop contributions 
from the scalar-neutralino exchange Feynman diagrams, for instance, involve the invariant 
misalignment parameters, \B A v\ and |S A A|. By contrast, the lepton number conserving 
observables can also depend on the scalar products ■uj = J2a''^% v ■ B — J2a''^aBa, and 
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To avoid the excessive effort involved in developing the basis independent formalism, we 
have pursued the analysis by making a fixed choice of the La field basis, while distinguishing 
between the free and constrained parameters phenomenologically. Although the choice f « = 
is perfectly admissible, this biases the exploration of the parameter space. Upon working 
in the basis choice Vi ^ 0, the constraints from the tree and loop level contributions to the 
neutrino mass matrix can be implemented in several ways. Since the tree contribution is 
proportional to X'^ = \ — Xvi/vd, we choose to regard Aj as free parameters and assign by 
hand Vi consistently with the restricted range of variation for A^. The basic relation between 
the gauge boson mass and the electroweak breaking mass scale, t>^ = + + vf = + v^, 
is then implemented by using the parameterization of the Higgs-sneutrinos VEVs 

Vu V Xj - X'i 

= 7 o = T; o TT^TT, Vi = GiVd, [ai = — r — , (11.28) 

tan/3 (l + (72 + tan2/3)2 A 

where we retain the usual definition for the ratio of VEVs, tan/3 = Vu/vd. At the one- loop 
level, the contributions to the neutrino mass matrix from the neutralino-slepton exchange 
Feynman diagram with a double mass insertion, as given by Fig. l(-B) and Eq. (111.54) 
of our previous work [4], bounds the misalignment parameter \B f\v\, or equivalently the 
effective parameters r]i= ^ — ^ — With the restrictions from the neutrino masses 

set on A^ and 77^, the parameters Vi and A^. are explicitly determined in terms of Aj. It is 
important to note that once the equations of motion are satisfied, the conditions A^ oc va 
and Ba = A^^^Xa oc Va entail the condition on the squared mass matrix itt-^j^^vb oc va- 
It then follows that fixing A^ and r]i still leaves the freedom of choosing the slepton mass 
parameters mf . Because of the strongly suppressed values A' << 1 and r}i « 1 imposed 
by the neutrino masses, the precise values assigned to A^ and rji have little effect on the final 
results. 

We shall develop the study of the regular vacuum solution of the scalar sector in terms 
of the neutral scalar field VEVs v^, f„, x and A necessary condition for stability is 
obtained by testing for the absence of saddle points of the scalar potential along the charged 
and neutral boson field directions in the field space. This condition is equivalent to the 
requirement that the squared mass matrices for the charged and neutral scalar bosons are free 
from negative squared mass tachyonic eigenvalues. To simplify the discussion we shall assign 
finite values for the lepton number violating parameters (coupling constants and VEVs) one 
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at a time for each lepton number flavor, so that finite couphng constants Aj, A^,, ^tt.^^^ 
and finite VEVs, Vj, are assigned for a fixed flavor index i. We thus neglect the quadratic 
terms of the form AjAj, ViVj, XiVj, ■ ■ ■, with i ^ j. The single dominant coupling constant 
assumption is not too restrictive as long as the lepton number violating contributions are 
small relative to those of the lepton number conserving contributions. To assess the validity 
of this approximation, one could envision using the variant prescription selecting some linear 
combinations of parameters in the SU (4) group of the La fields. 

The structure of the scalar mass matrices depends on the way in which one implements 
the equations of motion. Since the scalar potential of the NMSSM is given by a polynomial 
in the VEVs of the neutral fields, it is impossible, in general, to obtain analytic formulae for 
the said VEVs. Instead, we follow the practical procedure in which one eliminates the soft 
mass parameters rn^Hui ''^^h^i minimization conditions for the VEVs of electrically 

neutral fields Vu, Vd, and x, respectively. For the sneutrino VEVs, Vi, however, we need to 
specify beforehand our choice for the independent mass parameters, because of the presence 
of the off-diagonal and diagonal type sleptons mass parameters for the combination (Ha, !/«), 
namely and m| The freedom in solving the field equations (or equations of motion) 

for the sneutrinos Ui could be used to eliminate the mass parameters ? , which would then 
leave us with m| as free parameters. We call this option in the following as our prescription 
I. Alternatively, we could eliminate the mass parameters m| , and hence use ^ free 

parameters. We call this option in the following as prescription II. These prescriptions 
only differ in the way one treats the input data for the soft masses. While prescription II 
is perhaps more natural, since ni\, j are lepton number violating parameters on the same 
footing as the soft parameters Ay^ , this has the drawback of introducing inverse powers of the 
parameters Vi, which are expected to be small. In any case the results in prescription II are 
readily obtained from those in prescription I by substituting the expression for mf implied 
by the equations of motion for z>j, which we explicitly provide in Eq. (A. 9) of Appendix A. 
For definiteness, we quote below the formulae obtained with prescription I. 

C. Mass Matrices of Scalar Bosons 

The lepton number violating term proportional to Aj in the superpotential (A. 2) pro- 
duces mixing between the charged Higgs bosons and scalar leptons [4]. Similarly, there will 
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be mixing between the charginos and the charged leptons. These mixings can be studied 
through the appropriate mass matrices. Here we discuss the mass matrices of charged Higgs 
bosons/charged scalar leptons, and neutral Higgs bosons/ neutral scalar leptons, respec- 
tively. The mass matrices for charginos/charged leptons and neutralinos/neutral leptons are 
discussed in the next subsection. The field basis for the charged Higgs-slepton scalar modes 
is denoted by the column vectors 

= {H, \ e7, el), = {H^, g„ ef), (11.29) 

with the mass term in the Lagrangian given by 

-Lmass = ^Ih^lh^ch + H. C. (11.30) 

After a straightforward calculation, we obtain the mass squared matrix for the charged 
scalars whose elements can be written as 

,2 



,'^1 / 4 2,2 2 2 2 22 22, a\2 2 , 4\2 2\ 
+ TTli^'^L- + 9lVd -92Vd -9lVu -92Vu + "^X^Vu + "^X^x ) 



- ' 2 , „2 



+ 7^(fi'i + fl'2) + XdjiXdkiVjVk, 



2 

f2 _ 92'^dVu x2,. , \„/ A I \\ 2 



^H^Hu = 2 ^ ^"^^^ ^X{^>^ + I^X) - XXiVuVi - jJLdu, 

2 

M2 \ \ 2 I 92^dVi , , 

H.e* = XXiX H XdAjXsijVAVB, 

' Z 

^Hael = -XdAiXAVuX + A^^^XdAiVA, 

MliHu = ^{-2i4^ + 9lvdVu-2X^VdVu + 2AxXx + 2KXx^) 

+ — (-/^L - '^XXiVdVu + xXi{A~^^ + Kx)) + ' 



^ije* = -/^L - ^^K'i^dVu + 2;Ai(Ax. + kx) + ^(^2 - 2A-), 
1 

^lic* = + ^[(S'l - 9l)iyd^ - vj) + (S'l + 5'2)^i^ + 4A-a;^]5ij + XAikXsjkVAVB, 



11.31) 

H.32) 

H.33) 
H.34) 

11.35) 

H.36) 
H.37) 
n.38) 
H.39) 

11.40) 



M~.~c = XaVuXaijX - Aj^ijXAijVA, 

M|«gc = m|, . . + ^9i{-Vd^ - Vi + Vu)5ij + XabxXab'jVbVb'- 

For completeness, we have quoted the results for general bases of the slepton fields. The 
above formulae agree with Ref. [34] for the MSSM with R parity violation. 
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We now consider the neutral spin-0 Higgs-sneutrino fields whose real and imaginary parts 
carry positive and negative CP quantum numbers. The field basis for the CP-even and CP- 
odd Higgs-sneutrino scalar modes are defined by the Lagrangian mass terms as 

= ^UmI,^^,, + ^lMl,^%,j, (11.41) 
where the neutral spin-0 fields have the decomposition 

^neut,i = -^{^s,i + i%,i), ^neut,i^{H% S, i>j). (11.42) 

The mass squared matrix for CP-even scalars in the prescription I, using the equations of 
motion for i>i to eliminate the mass parameter i^^^i. with m| as free parameters, is given 

by 



1 V OCX) ' ~ 

— {^G\vi + Vux{Ax\ + k\x)\ - ^^{A-^.\ + k\x) 
+^{ml + 2Gl{vd' - vj) + ~XKvu' + x')) + 

Vd ' Vd^ 

^s,du = -4:GlvdVu + 2X'^VdVu - AxXx - kXx"^ + 2XXiVuVi, 
Mlds = -AxXvu + 2X{Xvd - Kv^)x + 2XX,xv,, 

Mldu, = -Kx{A-^~Xi + KXix) - (m|^ - 2Gl{vd^ + v^^) + Xf{vu^ + x^))vi 



Vd 



-2GlvA 



s,uu 



= —[iG'ivu^ + Vdx{AxX + kXx) + x{A-^_Xi + KXix)vi], 

Vu 

= -AxXvd + 2{-{KXvd) + (A^ + X'f)v^)x + {-A-^.Xi - 2/^A^x)^;^ 
= 2XXiVdVu - x{A~^.Xi + KXix) + 2{-2Gl + Xf)vuVi, 



^s,Si>i = -A~^,XiVu + 2XXiVdX - 2KXiVuX + 2X'fxVi, 
Ml,^,^ = ml + 2Gl{vd' - v^') + ^{v^' + x') + QG^, 

and for the CP-odd modes by 



11.43) 
11.44) 
11.45) 

11.46) 
11.47) 
11.48) 
11.49) 

11.50) 
11.51) 
11.52) 



'ux{AxX + kXx) Vux{A^,Xi + KXix)vi 



2 

^ ^d ^ 



{ml + 2Gl{vd' - vj) + X]{v^ + x^)W 2G\v,' 

+ — 5 + — (n.53) 
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m;^^„ = x(AaA + «Ax), (11.54) 
KdS = Vu(^aA-2«:Ax), (11.55) 



(11.56) 



Vd 

2 _ VdXjAxX + KXx) ^ x{A~^.Xi + KXiX)Vi 

K,us = M^x^ - 2«Ax) + (A'^JXi - 2KXix)vi, (11.58) 
Ml^,^ = x{A-^~X, + kX,x), (11.59) 

,.2 X , AxXvdVu , ^;„(A3^ A^ + 4/^A^x)^;^ 

-^p55 = (4«;AT;dT;„H h 3^««;a:) H , (11.60) 

^' X X 

Mlsu, = ^«(^A,A, - 2KXix), (11.61) 
= + 2G'^(^;/ - 1;^') + Xl{vJ + a:^) + 2G>,^ (11.62) 

The above formulae agree in the case of vanishing sneutrino VEVs, Vi — 0, with those 
obtained in our previous work [4]. 

The change of basis to the would-be Goldstone bosons in charged and neutral sectors is 
implemented through the transformation 

/ cos d — sin /3 \ 

. (11.64) 

Vsmp cosp / 

In order to project out the CP-odd modes and the massless Goldstone mode, G^, one 
simply needs to apply on the 2-dimensional Higgs bosons subspace the 2x2 matrix rotation 

{Hdi, Hui)^ {G^, A)'^ = TZp{Hfii^ Hui)^ 1 where the label / refers to the imaginary part 
of the fields, namely H^i = \/2^{Hd). The gauge basis for the physical CP-odd scalar fields 
is given by {A,'^s{S),'!^{i>i))^ — 'JZjslHfii, H^i, Sj)^ , where 7?.^ — diag{TZi3,l). In the basis 
with Vi 7^ 0, the zero eigenvalue Nambu-Goldstone neutral and electrically charged bosons, 
defined by the conditions M^ -^G^ — and M^f^ -^Gj' — 0, are given in our current choice 
of basis by the column vectors — {vd, —Vu, 0, Vi)'^ and G~ — {vd, —Vu, fj, 0)-^. The mass 
eigenstate fields are defined in terms of the mixing matrices which diagonalize the squared 
mass matrices as 

{^s,p;l)mass = [[/..pM^^t/J^ = {Ml^)diag]- (11.65) 
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Among the variety of observables of use in testing the scalar sector (partial decay rates 
Z ^ Hi + Hj, Z ^ Hi + Aj, Z ^ x° + x^) we mention, for later reference, the ratio 
of Z boson vertex couplings ZZHi to the Higgs bosons defined by ^zzHi = { ^gSu' Y — 
[{Uj)u cos P+{Uj),i sin ^f. 



D. Mass Matrices for Neutralinos and Charginos 

The mass matrix for the coupled system of neutralino and neutrino fields has been cal- 
culated in our previous paper [4] . Since the formulae for this matrix are rather complicated 
and were given in full form there, we do not reproduce these results here. We shall now 
concentrate on the mass matrix for the charginos/charged leptons. The mass term in the 
Lagrangian for the charginos receives contributions coming from the following sources: 

1. Contributions from gauge interactions: 



igV2mX'i,,(t>: + H.ci (11.66) 

where T" is the generator of the underlying gauge group. A" is the corresponding gaugino, 
and 0j, ipj are the components of the matter superfield. 

2. Contributions from the superpotential (in our sign convention): 



+ - 



+ H.C. 



_d(f)id(f)j 

where W is the superpotential, and 0j are the scalar components of a chiral superfield. 
3. Soft supersymmetry breaking gaugino masses: 

1 



(11.67) 



(11.68) 



where"2" here refers to the SU (2) l gauge group, and i are SU (2) l indices. Putting together 
all these contributions, we can write the mass term for the charginos/charged leptons as 



Ml 







(11.69) 



where we have chosen the basis 



X^^ = (-«A+,^+,e;^,/i^,r^), 



(11.70) 
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X 



with the mass matrix for the charginos/charged leptons given by 

/ M2 g2Vu \ 

g2Vd Xx (A^)iiWi (A^)22^^2 (A^'jaa^^a 
g2Vi Xix -{X'')iiVd 

g2V2 X2X -(A'=)22t'd 

\g2V3 X3X -(A'=)33^;<i/ 

Here we have assumed that the lepton Yukawa couphng and mass matrices, A^^, 
are diagonal 



M| = 



(11.71) 
(11.72) 



(11.73) 



-ML/vd, 



III. IMPLICATIONS ON SCALAR SECTOR AND VACUUM STABILITY 

Before discussing the scalar sector, we briefly discuss some numerical predictions for 
the light neutrino masses. The tree level contribution and the supposedly dominant one- 
loop contribution arising from neutralino-slepton exchange, give a neutrino mass matrix 

of approximate structure in the lepton generation space, m^^ij ~ aA'^A^- + brjirjj, [A^ = Aj — 
Xvi/vd, Tji = B'JB = Bi/B — Vi/vd] where the coefficients a, b are expected to have a smooth 
dependence on the parameters A, tan and x. The pair of non- vanishing mass eigenvalues 
of this matrix are expressed in terms of the invariant parameters by 

1 



[A'^ = EV, v' = j:vh |A'A^r = -E(A:77,-A;,7,)^ ^'■v = J:>^M- (nu) 

i i i,j i 

For clarity, we quote the explicit expression for the single non- vanishing neutrino Major ana 
mass eigenvalue present at the tree level [4] 

ax'V^{M2/M, + gl/9l){K - ViX/vdf 



m^3 = 



y„2(M2/Mi + gl/gl) ( 4kAxV tan (3 + (p/ tan (3 - Xv^Y ] - M2A%V 



tan2/3 

[a = 2{kx - — ), p = -Y^XaVa, a* = -Xx, K = ^^]- (III.2) 

X ^ \/ 1 

We have attempted to extract the values of the coefficients a and b from numerical 
calculations of the tree level formula, Eq. (III.2) , and of the one-loop amplitude in Eq. (III.54) 
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of our previous work [4] . Using the reference set of parameters given in the caption of Fig. 1 
with tan/3 = 2, gives us the following order of magnitude estimates: a ^ 10 GeV, b ^ 10~^ 
GeV. For these predictions to be compatible with the limit on the neutrino mass, \m,^\ < 1 
eV, we need to impose < 10^^, rji < 10~^. We have not attempted a global fit to the 
neutrino experimental data but note that this could be done by following a similar treatment 
as that used for the MSSM with bilinear R parity violation [35]. 

A. Numerical Study of Scalar Sector 

In order to have a rehable semi-quantitative description of the NMSSM, we have added 
to the tree level scalar potential the lepton number conserving one-loop contributions from 
the top quark and stop pair of squarks. The corrections to the neutral scalar mass matrix 
are obtained by evaluating Eq. (11.25) using a standard procedure [32], so we refrain from 
quoting the explicit formulas. The results of our numerical calculations for the mass spectra 
of the scalar and sfermion modes are shown in Figs. 1, 2 and 3 as single parameter plots 
displaying the dependence on the lepton number violating couplings Aj. These results were 
obtained by considering the reference set of natural input parameters used in [32], which 
are detailed in the caption of Fig. 1. Although our analytical formulae include the trilinear 
lepton number violating interactions with coupling constants Xijk, X[ji^, in our numerical 
calculations we have retained only the lepton number violating couplings A^. 

In Fig. 1 we present the results of our numerical calculations for the masses of the lightest 
neutral and charged scalar modes, and neutralino and chargino modes as a function of the 
lepton number violating coupling, which we take to vary within the range Aj/A G [0, 1]. 
At small tan/3 = 0(1), the modifications induced by Aj are essentially quantitative. As 
the coupling increases, the lightest CP-even scalar mass is enhanced and the CP-odd and 
charged scalar masses are reduced. This effect of enhancement and reduction of masses are 
clearly seen at large values of r = x/v. These results are also sensitive to the value of the 
VEV ratio tan/3 = Vu/v^- While the constraints on the scalar masses restrict tan/3 to a 
small natural range in the NMSSM, this feature is even more pronounced for large values of 
the lepton number violating interactions Aj ~ A. Already for the value tan/i = 8, the mass 
spectra undergo a qualitative change with large negative contributions to both the neutral 
and charged scalar bosons causing the vacuum instability signalled by the occurrence of 
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tachyons at Aj/A > 0.5. This is explained by inspection of Eq. (11.28) which shows how 
increasing Vi/vd proportionately to Aj/A has the same effect as an increase of tan/3. The 
allowed intervals for tan further shrink with increasing x. We also present in Fig. 1 the 
masses of the lowest lying neutralino and chargino states. The results are characterized by 
a slow growth of the masses with increasing Aj/A. 

The coupling of Z to the Higgs, izzHn is found to vary in the intervals (0.84 
0.63), (0.99 0.91), (0.96 0.81), (0.99 0.94) as A^/A increases in the the range 
(0.1 1.0) for the mass spectra displayed in panels (a), (6), (c), {d) of Fig. 1. For 
the sake of a qualitative comparison with experiment, we note that the LEP data indi- 
cates a correlation between the lightest Higgs boson vector coupling and mass such that 
izzH^ < (0.02 - 1.0) for Mh, = (50 - 120) GeV. 

In order to understand the dependence of the scalar mass spectrum on the supersymmetry 
breaking interactions, we display in Fig. 2 the results obtained for the case A^, — 0. This 
choice is made only for illustrative purposes since it violates the alignment condition for the 
parameters rji. As compared to our initial choice, A^^ 7^ 0, in Fig. 1, we see that switching 
off the supersymmetry breaking trilinear lepton number violating interactions results in an 
enhanced effect at small tan (3 ^ 2 but a reduced effect at tan/3 ~ 8. Indeed, the vacuum 
instability at tan/3 = 8 sets in at larger values of Aj/A. 

In Fig. 3 we display one-parameter plots for the masses of the lowest lying neutral CP- 
even and CP-odd Higgs bosons as a function of tan /3 (fixed x) and as a function of x (fixed 
tan/3) for the discrete set of values for Aj/A = 0.1, 0.5, 1. From these results we see that 
the maximal values of tan/3 allowed by vacuum stability lie inside the range 5 — 10, and 
that this interval gets further reduced at larger values of the VEV of the singlet field x. 
As a related observation, we note that the allowed interval for r = x/v e [1, 6] shrinks to 
re [1, 2] at tan/3 = 4 for large enough values of Aj. Furthermore, we observe that for the 
fixed values of r = 2.03, 4.06 in the panels (a), (6), corresponding to xV2 — 500 GeV and 
1 TeV, respectively, our results at small Aj/A ~ 0.1 join smoothly with those at Aj = 0, as 
can be seen in Fig. 1 (panels (a) and (c)) of Ref. [32]. 
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FIG. 1: The masses of lowest lying CP-even and CP-odd neutral scalars, the charged scalars, and 
the neutralino and chargino states are plotted as a function of X-JX. We assume the single lepton 
number violating dominance hypothesis with a fixed generation Aj 7^ considered one at a time. 
The curves for Hj=i, Aj=i, Ci=i, x[Li SL^d Xi=i drawn in full, dotted, dashed, dot-dash and 
double-dot dashed lines, respectively, as illustrated in the legend. The upper panels (a) and (6) 
correspond to tan /3 = 2 whereas the lower panels (c) and (d) to tan /? = 8. The left hand panels 
(a) and (c), and the right hand panels (b) and (d) refer to the parameter r = x/v = 2.03 and 
4.03, respectively. Our choice of the parameters, as given in Rcf. [32] (Fig. 1), is described by the 
values: A = 0.5, k = +0.5, Ax = 500 GeV, = +250 GeV, and A^ = 1 TeV, thq = rufj, = 
1 TeV, Q = 300 GeV, which enter the tree and one-loop contributions to the scalar potential. 
In the prescription I in which we work, we set the sneutrino mass parameter at m| = 500 GeV. 
The neutralino-neutrino and chargino-lepton mass matrices are evaluated by assuming the relation 
between the gaugino mass parameters Mi = |i^M2, with /ci = |, ^2 = Ij while using the numerical 
values M2 = 400 GeV and Mi = 199.2 GeV. The values of parameters Vi and A^, are set through 
the alignment conditions A- 0, r/j 0, wl^ie A^ = Aj — r]i = ~ = [-^i = 

Bi-B^. Bi = A~~Xi, B = AxX]. 
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FIG. 2: The masses of the lightest CP-even and CP-odd neutral scalars are plotted as a function 
of Aj/A, for fixed Icpton generation, under same conditions as in Fig. 1. The single Icpton number 
violating dominance hypothesis is assumed such that finite coupling constants Aj of fixed generation 
are considered one at a time. We use the same set of NMSSM parameters as in Fig. 1, the 
only change being the vanishing trilinear supersymmetry breaking coupling, =0. We have 
also displayed the lightest neutralino and chargino masses although these are independent of the 
supersymmetry breaking parameter A^,, and are given by the same curves as in Fig. 1. 
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FIG. 3: The masses of the lightest CP-even and CP-odd neutral scalars plotted as a function of 
tan/3 at two fixed values of r = x/v = 2.03 and 4.06 (upper panels (a) and (5)), corresponding to 
xV2 = SOOGeV and 1 TeV, respectively, and as a function of r = x/v at the two fixed values of 
tan/3 = 2 and 4 (lower panels (c) and (c?)). We use the same set of parameters as in Fig. 1 and 
assume the single lepton number violating dominance hypothesis with a finite Aj of fixed generation. 
The three curves for the modes Hi=i (bold) and Ai=i (light) are drawn in full, dotted and dashed 
lines corresponding to the values Aj/A = 0.1, 0.5, 1. 
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IV. CONSTRAINTS FROM UNBOUNDED FROM BELOW DIRECTIONS AND 
CHARGE AND COLOR BREAKING MINIMA 



A. General Considerations 

The stability of the regular vacuum with respect to the UFB directions and the distant 
CCB minima in the scalar field space are expected to give useful constraints on the NMSSM 
parameters. We briefly recall here the main features of the renormalization group approach 
to this problem which was extensively discussed in the context of MSSM [22]. This approach 
is particularly useful in the context of a gravity mediation of supersymmetry breaking in 
grand unified theories (supergravity GUT), where the unknown independent parameters 
consist of a small set of parameters and one can develop semi-analytic methods. The idea 
is to select suitable directions in the scalar field space along which the contributions to 
the scalar potential are dominated by the supersymmetry breaking terms while those from 
F— and D— terms are smaller and smoothly varying. In practice, the field directions are 
parameterized in terms of some single real variable, w, along which the scalar potential is 
given by a polynomial in w which receives, for a suitable choice of the parameter phases, 
negative contributions of 0{w'^) or 0{w^) from the supersymmetry breaking trilinear scalar 
couplings and mass terms, and positive contributions of ©(w^) or higher order from the 
supersymmetry terms. The higher order positive terms lift the potential upwards at large 
w. Even when absent at the tree level, higher order contributions to the potential, which lift 
it, always arise at the one-loop order. Requiring the scalar potential along these directions 
to stay above the regular minimum yields algebraic conditions on the coefficients of the 
polynomial. Considering, for illustration, the case of the field direction yielding the quartic 
order potential, V{w) — Aw^ + Bw'^ + C, where the coefficients A, B, C are known algebraic 
functions of the various parameters, we see that this is an UFB direction ii A = 0. At A > 0, 
this direction develops a minimum if i? < 0, as seen by evaluating the extremum with respect 
to w^. 



The regular electroweak symmetry breaking vacuum is stable against decay to the new vac- 
uum solution as long as V{wmin) > Vmin, where Vmin denotes the value at the regular 
minimum. An improved quantum stability condition can be deduced without a detailed 




B 




(IV.l) 
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knowledge of the one-loop radiative corrections to the potential upon invoking the logarith- 
mic dependence on the masses of states scaled by the running momentum scale Q. One needs 
only to assume that the above inequality involving the tree level scalar potential still holds 
but with the parameters replaced by renormalization group improved, momentum scale Q 
dependent, running parameters, 

Vufb{w; Q = Q) > Vmin{wo; Q = Ms), (IV.2) 

where the arguments w and Q along the field direction are set in the grand unification 
framework as w G [rriyyjMx] and Q = Max {g2W, \tw, Ms), and those for the regular 
minimum are set as wq = [vu, Vd] and Q = Ms, where Vu and Vd denote the the usual Higgs 
VEVs, and Ms the effective supersymmetry breaking mass scale. Extensive analyses of 
the vacuum stabihty constraints have been developed for the MSSM [22]. For the MSSM 
with broken R-parity, Abel and Savoy [23] inferred lower bounds on the trilinear R-parity 
violating couplings Xijk, X'^jf., and Hirsch et al., [24] discussed the implications on the bilinear 
R-parity violating interactions. 

Although a similar programme for the NMSSM with lepton number violation appears to 
be well motivated, its implementation is substantially complicated by the need to consider 
two-dimensional field directions which include the singlet field S — x as an independent vari- 
able, since there are no obvious correlations between the contributions from the electroweak 
singlet and non-singlet fields. Another complication stems from the fact that the renormal- 
ization group formalism is more complicated. This can be seen from the absence of studies 
of the dangerous field directions for the NMSSM beyond the preliminary works [36-38]. The 
progress achieved through the quasi fixed point solution for the running parameters [39, 40] 
or the numerical studies of physical constraints [41, 42, 44] are not of direct help to us in the 
present work. Since our main focus is on the lepton number violation, rather than solving 
the full-fledged problem, we shall follow a simple phenomenological approach, which we now 
describe. We assign natural values for the lepton number conserving Yukawa couplings and 
soft parameters which respect the regular minimum stability, and then examine the effect 
of increasing the lepton number violating parameters Aj, Bi, Vi, consistently with the con- 
straints from the neutrino mass matrix. This approach is similar in spirit to that followed 
in [24]. 

The field directions involve suitable subsets of the electrically neutral and charged scalar 
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states and the squark states. We continue using the single lepton flavor dominance hypoth- 
esis in which the nonvanishing parameters, Aj, A^,, Vi are finite only for fixed generation 
labels i, j, • " "• Useful intermediate formulas for the scalar potential of the electrically neu- 
tral fields and for the classical equations of motion of the electrically neutral and charged 
fields arc given in Eq. (A. 8), Eq. (A. 9) and Eq. (A. 11) of the Appendix A. 

To start with, we need the potential at the regular minimum. With our prescription of 
eliminating the dependence on the soft mass parameters iriff^, rUff^, nig, and by 
using the minimization equations with respect to v^, Vd, x, and Vi, the value of the scalar 
potential at the minimum is given by 



2\2 



VMIN(Vy,,Vd,X,Vi) = -G^(Vd -Vu ) 



GKXVdVuX +3X (Vu X +Vd (Vu + X )) 



+ x{—3AxXvdVu + x{3X^Vu^ — A^x + 3k'^x'^)) 

Vux{Ai,Xi 2hiXix) - 2XXiVd{vu + x^) 
+ \2Gl{-Vd^-rVu')-Xl{Vu'-rX^) 



V? - G\vt 



(IV.3) 



The above formula is seen to be independent of the slepton-Higgs boson soft mass parame- 
ters, so that it holds in the same form in the prescription where one eliminates m| rather 
than . We also note that the quartic and quadratic dependence on the VEVs coming 

from the gauge interactions combine into a term of same form as the minimum value of the 
scalar potential in the MSSM with bilinear i?-parity violation 



'(yMIN)MSSM = -G\{vl Vd - Vu) . 



(IV.4) 



B. Unbounded Prom Below Directions 

The Higgs-slepton field directions of interest are those which minimize the positive con- 
tributions from F— and D— terms and maximize the negative contributions from soft masses 
and trilinear couplings. There are three main unbounded from below directions defined by 



UFB - 1 
UFB-2 
UFB -3 



Hd = H^^O, S^O; 

Hd ^0, Hu^O, i>i^O, 0; 

Hd^O,H^^O, i>,^0,e, = e5^0, S^O, [i^j]. 



(IV.5) 
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The singlet field dependence of the potential makes an analytic study intractable. For 
instance, minimizing the potential with respect to S introduces a non-polynomial dependence 
with respect to the variables describing the electroweak non-singlet field directions. 

We now discuss the three UFB directions in detail. Along the UFB — 1 direction described 
by the two variables \Hd\ = \Hu\ = w and S = x, the potential is given by 

VuFB-i{w,x) = Aw"^ + B{x)w'^ + C{x), 

A = A^, B{x) = ml^ + m]j^ - 2Ax\x + 2i4^ + x\2X{X - k) + A^), 
C{x) = K^x'^ - \a^kx^ + mlx^. (IV.6) 

We have included the dependence on //^^ although this parameter is expected to be absent, 
as already discussed, in the minimal version of the model. For A 7^ 0, the would-be UFB 
direction is lifted at large and features a minimum at Wmin unless B{x) > It is 
useful to note that B{x) is a quadratic form in x with an extremum at Xmin defined by 
dB / dx\x=x^in = =^ Xmin = ^A'^/[2A(A — k) + Xf]. Since this is a minimum provided 
that 2X{X — K) + Xf > 0, one can express the condition for the absence of the UFB-1 direction 
by the appoximate bound on the soft parameters 

B{x) > B{xrmn) = + - ^^^^^^'^^^ + > 0- (IV.7) 

It is interesting to compare with the corresponding bound in the MSSM, B{xmin) ''^'h^ + 
m|^^+2(/i^— (A^/i)^). To obtain the improved constraint, one should determine the position of 
the potential minimum, Xmin and Wmm, and require the condition VuFB-i{xmin, Wmin] Q) > 
Vmin{x = xo, Vu^v sin (5,Vd^v cos Q = Ms). 

The UFB-2 field direction can be conveniently described by the parameterization Vu — 
w, Vd — wcosO, Vi — wsinO, x, [9 e [0, 27r]], which is designed to cancel the D— terms. 
The scalar potential is given by the quartic order polynomial in w 

VuFB-2iw,x,0) = A{x,0)w^ + B{x,0)w'' + C{x,0), (IV.8) 

where 

A{x,e) = (A cos(^) + Xi sm{9))\ (IV.9) 
B{x, 0) = X^x^ + A- x^ + -2x {AxX + nXx) cos{0) 

+ (A^ x"^ + cos^(6l) -2x {A~^.Xi + KXix) sm(0) 
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+ (\^x2 + m|.)sin2(e) + AAiX^ sin(2e), (IV.IO) 
C{x,e) = ''^^f' +K'x' + x'ml (IV.ll) 

Since the UFB-2 direction coincides with the UFB-1 direction at ^ = 0, to determine 
whether one avoids a distant minimum along w it is only necessary to test the condition 
that the coefficient B{x, ^) < at some finite 6. Inspection of Eq. (IV.IO) for B{x, 9) shows 
that the condition B{x,6 — > is sensitive to the signs of parameters rn^j^, m|^^ and Ax, 
while the condition B{x, ^) > at finite 9 is sensitive to the signs of A~^, and m| . To achieve 
B{x, 0) > 0, it is more favorable to restrict to the choice Ax < 0. We now attempt to assign 
by hand typical values to the relevant free coupling and mass parameters, and determine 
numerically whether the lepton number violating interactions can drive 9) to negative 
values if it started from a positive value at = 0. The dependence on 9 is displayed in 
Fig. 4 for typical values of the input parameters which are specified in the caption of that 
figure. We see from these results that for trihnear couphngs A~^, of same negative sign as 
^4^, the Aj interactions give significant positive sign contributions to B{x, 9) which avoid 
the occurrence of dangerous minima along w at finite 9. By contrast, choosing triscalar 
couplings of opposite sign, Ay^_ > 0, always drives B{x,9) to negative values at finite 9. We 
note that a dependence on 9 of similar type is found with various other choices of x and that 
changing the sign of Aj leaves the potential unchanged up to the replacement, 9 ^ tt — 9. 
We, thus, conclude that the lepton number violating couplings, Aj and A^, , have the ability 
to remove or induce the dangerous minima at finite 9 depending on whether A^, is negative 
or positive, namely, of same or opposite sign to Ax- We note that the case with Ax and A^, 
of opposite signs clashes with the universality of supersymmetry breaking and is unlikeky 
to occur in the context of grand unified theories. 

The UFB-3 direction can be conveniently described by the parameterization = 

aw'^, Cj = = (^w, Vi = aw{l + w"^)^ [i ^ j], designed to cancel the D— terms, while 
using cr = 4|, in order to satisfy the F— term fiatness condition Wh^ = 0. The potential 

3 

along this direction is given by 

VuFB-3{w,x) = -\a^kx^ + k'^x^ + \] + \] + (-2 A-^~Xi Vl + w^x 
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B(x =100, theta), l_i/l=(0.1, 0.5, 1.) 




FIG. 4: The coefficient B{x,6) for the UFB-2 direction is plotted (in GeV^ units) as a function 
of 9 for the relevant input parameters set as, A = 0.7, k = 0.3, Ax = —100 GeV, m|^^ = m| . = 
+100 GeV^, mjf^ = -100 GeV^, x = 100 GcV, and a discrete set of choices for Aj and A^^. 
The three curves for A^, = —100 GeV in the three cases Aj/A = (0.1, 0.5, 1.) are drawn in 
hght fines with dashes of increasing lengths. The three curves for Aj/A = 0.5 in the three cases 
A^ = (100, 200, 300) GeV are given in thick lines with dashes of increasing lengths. 
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-2 K Xi VTTuPx'^ a^) + w'^ (7^ (2 A- x'^ + m|. + m|e) 

(7^ (A^ + 2 A- + + m^J + m^x^ - 2A^^jXijja^'w'^{l + w^)^ 
+(7=^^/;^[2Aa,-,Aax + A^^.^.a + X^^^aw^l + w') + Aj^.(7(l + w% (IV.12) 

Inspection of the above potential indicates that the lepton number violating contributions 
from Aj or Xijk add positive terms to the potential that have the ability to lift the UFB-3 
direction. We restrict ourselves to an illustrative example by assigning the NMSSM parame- 
ters typical values resulting in a potential with a deep minimum at Aj = 0, and determining 
whether switching on Aj to finite values lifts this minimum. In Fig. 5 we display a repre- 
sentative case in which a deep potential well, produced by choosing negative squared mass 
values for rn\j^ and gets removed upon increasing Aj/A. We conclude that the lepton 
number violating interactions can be effective for lifting the UFB-3 field direction. 



C. Electric Charge and Color Breaking Minima 

The lepton number violating interactions may generate a minimum of the scalar po- 
tential along the field direction involving finite VEVs for the charged Higgs-slepton fields, 

—< Hj >, Ci —< Ci >. This so-called type II charge breaking minimum has been initially 
considered by Hirsch et al. [24] for the MSSM with bihnear R-parity violation. We pursue 
a corresponding analysis for the NMSSM with lepton number violation by solving the equa- 
tions of motion for the neutral fields in terms of the soft mass parameters iri'jj^, fnjj^, m| 
and m| (corresponding to prescription II), and substituting these into the two equations 
of motion for the charged fields, V- and e^. The equations are displayed in Eq. (A. 11) of 
Appendix A. Note that we can ignore the VEV of the field since — appears to be 
the only solution for Af = 0. The equations for v_ and Cj depend on the gauge and Yukawa 
couplings, the soft mass parameters ^ , and Vi. 

To start with, we consider the limit of small lepton number violating couplings obtained 
by expanding the equation of motion in powers of the small parameters Aj. Substituting the 
expansion up to leading order in Aj 
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FIG. 5: The potential energy density V{x,w) for the UFB-3 direction is plotted in GeV units 
as a function of w for the choice of input parameters, A = 0.7, k = 0.3, Ax = 500 GeV, = 
500 GeV, Ai = 500 GcV, rajr = mjr = m\ . = -200 GeV^ = ml = 100 GcV^, with 

' Aj ' rid L.l ^ ^ ' 

VEV parameters, x = 100 GeV, tan/3 = 2 and a = \x/\% = Xxv^jmr-, [m^ = 1.777 GeV] for the 
tau-lepton field case. The reduced range of variation of the variable w is explained by the large 
value assumed by the scaling factor, a ^ 30.6 x. The curves for Aj/A = (0.1, 0.5, 1.) are drawn 
with dashes of increasing lengths. 
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Vi = Xivf-^ , = ~XiV^y , a = ^ , (IV. 13) 

in the equations of motion for V- and e^, we find tliat the equations in leading order of Aj only 
admit the trivial solution v- — 0, Ci — 0. Thus, we reach the conclusion that as long as the 
couplings Aj are small compared to unity, the lepton number violating interactions cannot 
cause the emergence of charge breaking vacuum solution. This conclusion for NMSSM 
corresponds to that reached in Ref. [24] for the MSSM. 

To determine whether a non-trivial solution is favored at finite Aj, wc attempt to solve 
numerically the equations of motion for V- and Cj, and compare the minimum value Vmin 
with the regular minimum value Vmin at a discrete set of values of the parameters A^ and Vi. 
For fixed tan (3, we evaluate Vmin by using the formulas Vu — v'sm(3, va — v' cosfi, [v' — 
{v'^ — v^)^. Based on the argument that the singlet VEV x is not strongly affected by the 
lepton number violating interactions, we identify the value of x in Vmin with its value along 
the charge breaking field direction in Vmin- The values of tan j3 near unity arc critical, as non- 
trivial solutions exist only around tan/3 = 1 and disappear quickly at larger values. Since 
the task of determining the exclusion plot in the Aj, Vi plane is cumbersome, we use typical 
values for the relevant parameters, so as to determine to what extent a charge breaking 
solution at Ai = becomes disallowed by increasing Aj. We restrict ourselves to solutions 
with real values of v_ and e^. For the choice of input parameters A — 0.7, gj = 0.127, g'l = 
0.425, K = 0.3, Ax = 500 GeV, A, = 500 GeV, A^ = 250 GeV, ml . = 100 GeV^, 
with the fixed values of the VEVs x = 100 GeV and Vi = v/10, we find the real solution 
{v_,ei) = (-60.9,585.1) GeV having Vmin - Vmin ~ -8.26 x 10^ GeV^ for values of 
tan/3 = 1 and Aj/A = 0.1. However, we find no non-trivial solutions as we increase the 
couphng Aj/A e [0.5,1]. We thus conclude that the lepton number violating interactions 
have the ability to lift the charge breaking minima. 

Finally, we comment briefly on the issue of charge and color breaking minima in the 
NMSSM by focusing on the field directions [18] described hj Vu = Ui = = w and the 
singlet field VEV x, with all other fields vanishing. Assuming x to be frozen, for simplicity, 
one finds that the resulting potential, 

VccB = K^x^ + mW - Ia^kx^ + w^[3Xfw^ + 2A^A> + ml, + ml^ + mj,^ + 

6 ' ' A 

does not develop a deep minimum along for small Yukawa coupling constants A", 
provided one satisfies the conditions on the trihnear scalar matter couplings, A^ < 
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3(m^^ + m~. + + (A^a^)^)^- This approximate result suggests that the bounds on 
Af should become weaker upon increasing A,. The general field direction described by [22] 

\ui\, \Ui\, \vd\, \vj\) = (1,q;,/9, 7, 7l)w along with the singlet field VEV x, involves a 
more elaborate discussion which we shall not pursue here. 

V. CONCLUSIONS 

In the present work we have examined the effect of lepton number conserving and violating 
Yukawa couphngs of same size, Aj A, on NMSSM. One expects significant modifications for 
the scalar sector observables since the spontaneous electroweak gauge symmetry breaking 
is now linked to both the down type Higgs boson and sleptons. An important challenge 
was raised by the need to define a simple parameterization of the model consistent with 
the constraints on the light neutrinos without specifying in detail the underlying dynamics. 
Noting that the dominant contributions to the neutrino Major ana mass matrix are controlled 

~ ~ A- A 

by the effective alignment parameters A^ = Aj — Xvi/v^ and r]i — ~ ^' ^® proceeded by 
treating Aj as free parameters while determining the values of the sneutrino VEVs Vi and 
Aj^_ through the restrictive conditions, A^ ~ 0, T^j ~ 0. There is no unique prescription, and 
the one used may well single out a non-generic region of the NMSSM parameter space. To 
answer this objection, one could attempt building a supergravity unified model with a 
symmetry spontaneously broken in a hidden sector so as to check whether this respects an 
approximate dynamical alignment robust under the renormalization group scale evolution, 
following a similar analysis as that of Nilles and Polonsky [27] . 

In the first part of the present work we examined the impact of the lepton number 
violating interactions on the mass spectra of scalars. The vacuum stability constraints were 
found to restrict the VEV parameters, tan (3 — v^/vd and x, to narrower intervals than in 
the lepton number conserving case. This property is refiected in the fact that the bounds 
on Xi become stronger for larger tan (3 and x. There are certain analogies between our study 
and that developed for the MSSM with bihnear R-parity violation by Davidson et al. [14] 
in terms of the basis invariant parameter 5r, corresponding to J2i Vi basis choice, at 

fixed values of the CP-odd scalar masses. The conclusions in the latter work regarding the 
reduced range of variation of tan P and the size of the corrections to the scalar sector masses 
are qualitatively similar to ours. However, no meaningful comparison can be made because 
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of the different parameterizations. 

Tfie second part of the present work was devoted to a quahtative study of the vacuum 
stabihty contraints from the UFB field directions and the CCB minima. The discussion 
for this case, unhke that in the first part, does not bear directly on the implementation of 
the parameter alignment conditions. Based on illustrative examples covering a small part 
of the parameter space, we found that the lepton number violating interactions may have 
a positive impact on the regular vacuum stability provided the coupling constants Aj and 
A^, assume large enough values. Since the approach of selecting field directions based on 
the renormahzation group cannot be developed in the NMSSM by analytic means only, its 
advantage over a systematic numerical exploration covering the full field space is not clear. 
Nevertheless, our discussion indicates that pursuing the renormahzation group approach 
on more quantitative grounds is worthwhile. The recent progress in developing efficient 
numerical methods to search the global minimum of the scalar potential in multidimensional 
field spaces [45] could be useful for further studies along these lines. 

APPENDIX A: USEFUL FORMULAS AND CONVENTIONS 

The i?-parity conserving and i?-parity violating parts of the NMSSM superpotential can 
be written as 

Wrpc = X^i.HuQjU^ + Xji^HciQjDl + X'jf.HdLjE^ + XHdHuS - -S , (A.l) 
Wrpv = IXijkLiLjEl + KjkLiQjDl + ^XJ^.U^Dpi + ~\L,H^S. (A.2) 

The covariant four vector notation for the lepton and the down-type Higgs superfields, and 
for the Yukawa couplings, employed in this paper is 

Xa — (A,Ai), XAjk — {^"jk, Kjk), ^Ajk — i'^jk^ Kjk)- (-^-4) 

We use the following convention in writing the multiplication of chiral electroweak doublet 
superfields: 

LaH^ = LA-e-H^^ vaHI - caH^ = (i^X - H^H+) + {v,Hl - e,i/+), (A.5) 
H^Q = H^-e-Q^ H+D - H^U, HaQ = Ha-e-Q^ H^D - H^U, (A.6) 
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where e is the 2x2 antisymmetric matrix with eu — —621 = 1. The it!-parity conserving 
and i?-parity violating contributions to the trihnear part of the potential can be written as 

ft ~ ~ d d ~ ~ ~ ~ A K 



3 

V^-J^ = -\A^,k\jkUL,El - A^;,X^,MjDt - lA^kKjkU-Dpi - A-^^L^H^S + 4Ac7.) 



The relationship that we adopted between the superpotential and the effective Lagrangian 
uses the convention Leff = + H.c. = — + \Wijipii/jj + H.c. in contrast to 

certain authors which use the opposite sign convention, Leff = —{[W]f + H.c.) = — — 
jWijipiipj + H.c. The latter sign convention is that adopted, for instance, in the studies of 
the NMSSM by Miller et al, [31] and by Barger et al., [32]. Accounting for this fact, we 
obtain the following correspondence between our notations and that of the latter authors: 
A hs, Ax — > As, K ^ K, — >■ ~A^, A" ht, A^ At, a; — > s/\/2. We have compared 
our formulas for the scalar sector potential and mass matrices and for the neutralino sector 
and found complete agreement. 

The scalar potential for the electrically neutral fields, with the electrically charged and 
color non-singlet scalar fields set to zero, is given by 

Vp + VD + Vsoft = \XVd+~XiVi\'^\x\'^ + {\Xvu\'^ +\~XiVi\'^)\x\'^ 

+ \vuiXvd + XiVi) - Kx^p + - Ivdl"^ - \vi\'^y 

~ A K 

- AxXvdVuX - A-^XviVuX '^\x\^ + ii\yVAVu + 'rn\i'LVdvl + H. c. 



+ 

+ mfjjVu\ +rnjjj\Vd\ + mi,\Vi\ +ms\x\ . (A.8j 



The relevant neutral and charged field VEVs are denoted as, < H^ >= Vu, < H^ >~ 
Vd, < S >^ X, < Vi >= Vi, and < H^ >= V-, < >= e^, < >= e?, < qi >= qi, < 
ql >— qf, [qi — {ui, di), qf — {u^, c?^)]. Our choice of field basis, v+ — 0, obviates the 
need to consider the electrically charged direction < if+ >= The formulas determining 
the soft masses of the Higgs bosons and sleptons through the minimization equations of the 
neutral fields in prescription I (using m| as free parameters) for finite values of the neutral 
and charged fields are given by 

2 fi'l('2i2i2 2\ , 92 / 2 2,2 2\ \2 2 

rriR^ = ^{ei+v^+v^- w„) + — (-e^ - v_ + - v„) - A 

^^IaAW _^2^2 _ S^2 ^2 I I^XVdX^ 
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+{-2X~XiVd + 



2 9\ I 2 2 2 I 2\ I fi'2/' 2 2 2 , 2\ 



\2 2 , ^A^t'u^; 

-A v„ + 



X^x'^-r 



Vd 



4 

K A v,, x"^ 1 



{-A-^,XiVuX - nXiV^x^) Vi 



I I 9l f 2,2 2\fl'2 / 2 2 2\ I 2 , \2/ 2 , 2\\ 2 , 9l'^92 4 



m| = —vJ^ A^ — 2 ej v_ A Aj — e^^ A^ — A^ + 2 k A v „ — A^ 



+ Ai^K X - 2k^ x^ + (-2 A Aj + 2 K Aj v„ + 



Ak Xi Vu 



X 



)vi-XiVi , 



1 ~ 1 

^^dZ, = -2^2 t;- - A Ai (^;^ + x') + — 



XiVuX{Ay + KX) 



+ 



9l / 2 2 2 , 2\ , 92 / 2,2 2 , 2\ 2 \ 9l ^ 92 3 

-ji-ei -v_-va + vj + — (-e^ +v_-Va + vj - m-^,jvi — v? 



(A.9) 



The equations of motion for the La — {Hd, Lj) fields can be expressed in the SU (4) group 
covariant notation as 



^1,ab'"b = XaXVu{Ai^ + KX), 



[Ml 



L,AB ""LAB 



irif 



+ XaXbK + x") + 2GX{vl - vI)5ab\ 



(A.IO) 



The above structure of the neutral scalars squared mass matrix satisfies the important 
property that the alignment A^ oc va is satisfied if and only if A^ is an eigenvalue of the 
matrix M| This result generalizes that found in the MSSM with R-parity violation [10, 
46]. 

The minimization equations for the charged slepton fields in prescription II (using 
as free parameters) are given by 



el: = 2 ef gl + e'i - [e^ g'f + 2 m|. - g't {vJ + v^' + Vd' - vJ)] - 2 A^ A^e, Vd, 



Ci c/2 V- Vd + v_ Vu {g2 VdVu-2X VdVu + 2 AxX x + 2kXx ) 



= — 

Vd 



+2 Ci {m%^i, Vd + XXi Vd x'^) 



+ - 



Vd 



[Ci gl {-vJ + Vd') -2v_ {ml^-^^ + XXi {vu' + x'))) v, 



9l V-Vi', 



ej : = — -{eig2V^Vd + 2 {m]j^i, Vd + XXi Vd v^ - A~^.Xi VuX + XXiVdx'' - k Xi v^ x^)) 
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+ei (-2 A- Vu + gl {vJ - Vd' + Vu)) + 2 1;_ + A Aj x^) + gl v_ 



(A.11) 
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